Glassy trapping of manifolds in nonpotential random flows 
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We study the dynamics of polymers and elastic manifolds in non potential static random flows. 
We find that barriers are generated from combined effects of elasticity, disorder and thermal fluctu- 
ations. This leads to glassy trapping even in pure barrier-free divergenceless fiows v '^^Z" {(f) > 1). 
The physics is described by a new RG fixed point at finite temperature. We compute the anomalous 
roughness i? ~ and dynamical t L" exponents for directed and isotropic manifolds. 
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Elastic structures in presence of quenched disorder ex- 
hibit glassy behaviour such as pinning and very slow 
thermally activated motion of velocity v ^ e^^'^^'^^'^ . Un- 
derstanding these properties is crucial for numerous ex- 
perimental systems, e.g. vortex lattices in superconduc- 
tors, magnetic domains, solid friction and charge den- 
sity waves For £) < 4 the static system is usu- 
ally in a disorder dominated glass state described by a 
T — Q RG fixed point Q| . Barriers diverge at small force 
U{f) ~ /"^ and large scales §. If D < 2 thermal fluctu- 
ations become important again. Periodic structures e.g. 
still exhibit a glass phase in _D = 2, but with weaker log- 
arithmic barriers v ^ described by T > RG fixed 
points . Up to now, the study of pinning and glassy 
dynamics has focused on potential systems, typically an 
elastic manifold in a random potential Q . The existence 
of pinning and barriers in a non potential system was pro- 
posed recently in the context of vortex lattices driven over 
a disordered substrate Q . Driven systems even when all 
microscopic forces are formally the gradient of a poten- 
tial, are usually not so at a macroscopic level since energy 
is constantly pumped into the system. In ||^ it was found 
that despite the non potential convective term present in 
the equation of motion, some glassy properties survive. 
Related observations were made in (mostly mean field) 
spin models H,^. Since little is known about general 
properties of non potential systems with disorder, it is of 
great interest to study other examples (and beyond mean 
field). An outstanding question is whether the effective 
temperature, which tends to be generated in these sys- 
tems from the constant dissipation, does overwhelm the 
effect of static disorder. 

In this Letter, we study a model of a polymer in a 
non potential static random flow (a "randomly driven 
polymer" ) and its generalization to an elastic manifold. 
We work directly in finite dimension using dynamical 
RG. We study two cases: either the polymer is directed 
(each monomer sees a different flow) or it is isotropic (all 
monomers are in the same flow) see Fig. 0. Our main 
finding is that the physics is described by new RG fixed 
points with both T > and finite disorder, leading to 
anomalous roughness ( and glassy trapping by the flow 



(with sublinear v{f) /*). Compared to the single par- 
ticle studies [|lO| , new features emerge such as the crucial 
role of internal elasticity in generating barriers. While 
being consistent with the Hartree results {d oo) for 
long range (LR) disorder our present study yields 
new universal fixed points for short range (SR) disorder. 
It also generalizes the dynamics of self avoiding mani- 
folds to a quenched disorder situation, though self 
avoidance is perturbatively less relevant in most regimes 
studied here. For real polymers in realistic fiows ours 
is mostly a toy model which could be improved by in- 
cluding hydrodynamic forces. Nevertheless, some of the 
physics unveiled here will be present in more realistic 
models, such as the existence of preferred regions in the 
flows. Inhomogeneities in polymer distribution and non 
linearities in v{f) could be investigated experimentally. 





FIG. 1. elastic manifolds (polymers D = 1) in random 
flows: (a) directed polymer (b) isotropic chain. 

Let us illustrate how elasticity (or entropy) leads to 
dynamical generation of barriers in a divergenceless flow. 
It is well known that for a single particle convected in 
such a flow, the stationary measure at T > is spatially 
uniform and the drift velocity under a force simply v — f. 
Remarkably as soon as one considers two elastically cou- 
pled particles, preferred regions appear. Consider motion 
in the simple 2d fiow (vx = —y, Vy — x): 
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zi = c{z2 - zi) + i{uj + 5)zi + r]i 
Z2 = c{zi - Z2) + i(w - 6)z2 + ri2 
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Zk = Xk + iyk is the complex position of particle k, c 
the elastic coupling, < ry^yy* >= ATSij is the thermal 
noise. Without disorder 6 = the motion is just a 
collective rotation around the center (the matrix eigen- 
values are iuj and —2c + iuj). At T > there is free 
collective diffusion and V = f e.s before. As soon as 
S > the zero mode disappears and the two parti- 
cles converge towards the center (the eigenvalues are 
X± = -c + iLU ± \J(? - ?P-). This effect persists at T > 
since one finds < |zfcp 'lT{c^ + ^^)/(c(5^) = 1,2) 
and thus there is a genuine bound state despite the diver- 
genceless nature of the flow. The stationary conformation 
is twisted since < z^z^ >= 2T {c + iS) / S'^ . Extending (|l|) 
to a directed chain, one finds similar results (the decay 
towards the center at T = is exactly the problem of spin 
depolarisation Thus a directed polymer in a flow 

will be attracted to regions where elastic energy and dissi- 
pation are smaller (though the precise balance remains to 
be understood). Presumably polymers in realistic flows 
will be repulsed by high vorticity regions. Non linear 
extensions of the above model will thus show the genera- 
tion of barriers. Similar effects exist if all monomers see 
the same flow (isotropic chain). This is reminiscent of 
the effect discovered by Fauve and Thual mapping 
their Landau- Ginsburg function onto the complex posi- 
tion z{x) of an elastic chain in a 2d non potential non 
linear flow, one sees that (even at T = 0) it may converge 
to stable localized (rotating and twisted) conformations 
z{x,t) zo(a;)e*"*. 

We now turn to our full disordered model where these 
effects can be studied quantitatively. We consider a 
manifold of internal dimension D parametrized by a d- 
component fleld r"(x). The polymer corresponds to 
D = 1 (x labels the monomers), and a single particle 
to I? = 0. We study the Langevin dynamics: 

Vdtr^t = c^lrZt + [r^t , x] + (2) 

T] is the friction, c the elastic coefficient and the gaussian 
thermal noise satisfies (C^tCf'c) = 2r]TSai3S{t — t')S^ {x — 
x'). Angular brackets denote thermal averages and over- 
bars disorder averages. Fa[r,x] is a gaussian quenched 
random force field of correlations: 

Fa[r, x]Fi3[r' , x'] = Aapir - r')hx-x' (3) 

There are two main cases of interest. If the manifold 
is directed (e.g. a polymer oriented by an external field), 
then h^-x' = 5^{x~x'). If the manifold is isotropic, (e.g. 
a gaussian chain in a static flow), the force fleld does not 
depend on the internal coordinate F^ [r, x] = F^ [r] and 
hx-x' = 1. We consider a statistically rotationally invari- 
ant force field with both a potential (L) ( "electric" ) and 
a divergence-free (T) part ("magnetic"). We thus study 
correlators AQ,^(r) = /S.'^e^^'^ with, at small K: 

AQ/3 T^a-d/^ KgKp KgKp 

^ \9L 2 + 9T{0al3 7P2~)) i^) 



with Jj^ = J ^^^^ . We are mostly interested in SR cor- 
related forces a = d (force correlator decaying faster than 
r~'^). We also give results for LR correlated forces (force 
correlator decaying as with a < d). Note that SR 
random potentials Aap{K) ^ gi^KaK/s -\- gT(Saf3K^ — 
KaKp) (formally a > d) renormalize to the random force 
case, except in the potential case gT = 0, much studied 
previously ||J|,^ which is preserved by RG. 

We study the time translational invariant steady state. 
One defines the radius of gyration (or roughness) expo- 
nent {{rxt — I'x't)'^) ~ the single monomer diffu- 
sion exponent {{rxt — rxt'Y) ~ \t—t'\'^'^ and assumes scal- 
ing behaviour {{rxt — fmY) ~ \x\^'^B{t/x^) with C, = zv. 
The drift velocity under a small additional applied force 
/ in (|) is i; ~ {rxt)lt- We find i; ~ /"^ at small /, 
with > 1 and thus the polymer is trapped by the flow. 
Without disorder one has Co = (2-£')/2, vq = {2~D)/A 
for D < 2, and bounded fluctuations for D > 2. 
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FIG. 2. Regions in d, D plane where SR disorder is rel- 
evant for a directed manifold. Inset: Regions in a, d plane 
where SR or LR disorder is relevant for the directed polymer 
D = l. 

We start with simple Flory (i.e. dimensional) esti- 
mates. Balancing all terms in (0) (e.g. for directed man- 
ifolds it reads r/t ^ r/x^ ^ j.-d/2^-D/2^ yields: 

Cf = 24" = ^— ? , Cf° = 2i^r = (5) 
^ 2 + ^'^^^ ^ 2 + d ^ ' 

for directed and isotropic manifolds, respectively (with 
d ^ a for LR disorder). This suggests that the man- 
ifold will be stretched by disorder We now study 
(^,H) using the dynamical generating functional Z = 
J DrDfe-^°^'^''^^-'^"''^'''^^+'f^ with: 

So = J d^xdt tf^t{r^dt - cWDr^t - vTtr^t^r^^ (6) 
Sint = ~^ / {i^xt){if1't')^ai3{rxt - rx't')hx-x' 

^ Jx,x' ,t,t' 
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We used dynamical RG methods, both via a Wilson RG 
scheme (presented here) as in [||J^ and via a field theo- 
retical method using MOPE techniques as in p8[jTl| ] (de- 
tailed in The free propagator is B{x, t) ^~^{{rxt^ 
^00)^)0 a-iid the free response function is DR{x,t) — 
is invariant under the rescaling 
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(directed case) and 
(isotropic case) with a — d 
everywhere for SR forces. There are generically three in- 
dependent exponents z, C, (3) and from them one gets: 



i^^C/z , ^={z-0/{2-C + P) 



(7) 



using that / j'giD+z+a)i ^ ^ _^ y^{c-z)i yj^jgj- rescaling. 
Power counting at the Gaussian fixed point with no disor- 
der {z = 2, P — Q and C, — (2— D)/2), shows that disorder 
is relevant when (Fig. 2) d < d'^" for directed manifolds 
and d < cijf° for isotropic manifolds (resp. a < dc for LR 
disorder), with: 
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(8) 



Thus for directed polymers disorder is relevant for d < A 
and for isotropic polymers for d < 6. Power counting and 
symmetries show that the only relevant terms generated 
in perturbation theory are the one in (|^) . To lowest order 
the corrections to T, c and rj read: 
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2/2/i^if--'^(e-^'^fe'0) - e-^'^- ) (9) 



and [d - l)Sr] = {gL/_9T)v^T ( where Boo = J^.-^ is IR 



infinite for d < dc) |20|. To simplify the above expres- 
sions we have used that the free propagator satisfies the 
FDT relation 0{t)f^B{x,t) = TR{x,t). We note that a 
useful feature of the present model is that although the 
FDT relation is not obeyed by the exact Green's func- 
tions (since disorder is non potential) it is within each 
scale of the RG. The FDT violation will simply appear 
as a non trivial renormalization of T . The relevant cor- 
rections to disorder read: 
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K) (10) 

K-K'{B^,^+B^_^,) 



The small P expansion in (10) is well behaved and 
because of analyticity the only divergent terms gener- 
ated are proportional to Sap yielding 5gL — Sgx- One 
can thus take the limit P {K' = -K in (^). 
One checks explicitly from (|l^) that (i) starting with 
Aaf3{K) - gLKaK(3 + griSapK'^ - K^Kp) with both 



> and 5l > one generates SR random forces and 
that (ii) a divergence free random force field does not re- 
main so. The final RG equations have the form (specified 
to various cases below and in Table 1): 
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cdl 
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2-D- 2C + 13 + Argr 



^^eg^-iB- A)gLgT + Egl 
^^egj, + iA + E)gLgT - Bg^ 
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FIG. 3. Examples of RG flow diagrams, (a) directed 
manifolds with SR disorder. The physics is controlled by the 
fixed point /* at qt = Ql- (b) isotropic manifolds LR disor- 
der. The flow is along straight lines. There is a line of fixed 
points L and an apparent separatrix. C~Co changes sign upon 
increasing Ql/qt, suggesting progressive localization. 

(i) Directed manifolds, SR disorder, e = 2 — (2 — D)^- 
Then hy = 5^ {y) and there is no correction to the elastic 
coefficient 5c = (guaranteed by the statistical tilt sym- 
metry 1^ for SR and LR disorder) . Thus /3 = and there 
are only two independent exponents. The dimensionless 
couphng constant is g-^.L = CA<^gT,L/T('^+2)/2. The 
RG flow is depicted in Fig. 3. There is a globally attrac- 
tive isotropic fixed point /* at 5^ = 5t = ^li^ — A) 
(note that B > A) and the line = ^-j" is preserved. 
The potential line = is preserved with a flow to 
strong coupling. There is another apparent fixed point 
at gL = 0,gT = e/B, attractive along the line gL = 
(divergenceless flows) . However ioi D > this line is not 
in the physical domain since, as discussed above, a flnite 
gL term is generated and one flows to /*. At /* one gets 
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and V and cj) using (0). The expansion can be carried from 
any point of the line d* = 4/(2 — D*). Optimizing over 
D* as in ^ yields the estimates C = 0.625 ± 0.02, z = 
2.085 ± 0.02, (j) = 1.061 ± 0.015 (d=3) and C = 0.9 ± 0.1, 
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z = 2.2±0.1, <t) = 1.182±0.08 (d=2). For small e, C > Ci^ 
and Flory is thus likely to be a lower bound p3[ 



We thank L. Cugliandolo, J. Kurchan, T. Giamarchi 
and L. Schafer for discussions. 



LR disorder, e — 2 



(2-D)a 



(ii) directed manifolds, ±^±y. ^^.^^i^^,. ^ — ^ ^ 
with a < d. The novelty is Sg^^fi — since the LR part 
of the disorder is not renormalized |Tq]. The dimension- 
less coupling ijl^ Vl^t = CgL,T /T^"-^^^"^ will experience 
non trivial renormalization only because T renormalizes. 
Thus the ratio I^ty = -2^^ is preserved which leads to 
a line of fixed points. Thus = e/B and 

~^ ~9t- This yields the continuously varying ex- 



ponent z 



and C = Cf holds to all 



' ^QT (2+a)(d-l) 

orders in e in the LR case since neither the vertex nor c 
renormalize. As conjectured in is found identical to 
the Hartree approximation (formula (12) of ||^). 

(iii) isotropic manifolds (SR and LR disorder): Then 
e = D -f 2 - {a = d for SR disorder). The nov- 

elty is a renormalization of the elasticity of the manifold 
(it becomes stiffer 5c > 0) and thus a third non triv- 
ial exponent /3. The dimensionless constant is ^ — 
C 'il1'\~a ■ For SR disorder the RG flow is similar to 

the directed case. There is a fully attractive (isotropic 
disorder) fixed point /* at — '9l ^ 1/{B — E — A). 
As before, LR disorder is not renormalized, the flow lines 
are n\ = — = fc and there is a line of fixed points 
gT = e/ {B — kE), g^ — kgx parametrized by k. The gen- 
eral formula for the exponent is C = (2 — -D)/2 -I- e{AT — 
kAc)/2{B-kE-A), z = 2+ke{Arj-Ac)/{B-kE-A) and 
P — —keAc/ {B — kE — A), with the values given in Table 
1 (and set A: = 1 for SR). For SR disorder this yields for 
= 1, C = 0.5 + 0.13e, z = 2 + 0.04e, /3 = -O.OlSe. Ex- 
trapolations as in suggest ( = 0.85±0.1, z = 2.1±0.1, 
(3 = -0.025 ± 0.01 for polymers in d = 3. Note that 
for LR disorder the line of fixed points is not paral- 
lel to the axis any more (see Fig. 3), and thus there 
is no fixed point for k > kc — B/E. There is even a 
value k* = ^^kc < kc where C ~ Co changes sign. For 
k* < k < kc the localization effects dominate. Together 
with the apparent separatrix in the flow, this suggests 
two phases (a localized one and another one with a con- 
tinuously varying C). It would be interesting to check this 
scenario by non perturbative methods (such as in |^). 

Thus polymers (and manifolds) in static flows are de- 
scribed by non trivial RG fixed points at finite disorder. 
The main difference with conventional static (potential) 
glasses is that a temperature is generated which weakens 
disorder (and also renders the fixed point perturbatively 
accessible - the coupling constant being g jT^'^-^'^)/'^ !). 
The main similarity lies in the sublinear v{f) response 
reminiscent of the marginal glass in the statics Nu- 
merical simulations would be helpful to check the present 
results and further investigate the physics of this phase 
p^ . The competition between localization and driving 
found in some cases for isotropic manifolds (i.e correlated 
disorder) also deserves further studies. 
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A scenario to be checked is whether the hierarchy of wells 
within wells of a static glass is interrupted here at some 
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random level by the nonpotential forces (a mechanism to and if there are Levy walk type of jumps between trapped 

generate wells with algebraic distribution of release time) configurations. 
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Table I: Coefficients of the KG equation (|ff|) for cases studied 
here (one must set d = d* and D — D*). One denotes — 
/+°°dtB(f,t)-"/^ withB(l,t) = ^((4t)(2-°)/2e-i/(«) + 

r[D/2, l/(4t)]), with /d=6 ~ 1.7935, 



